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We show that fermion systems with random interactions lead to strong coupling between glassy order 
and fermionic correlations, which culminates in the implementation of Parisi replica permutation 
symmetry breaking (RPSB) in their zero temperature quantum field theories. Precursor effects, 
setting in below fermionic Almeida-Thouless lines, become stronger as the temperature decreases 
and play a crucial role for many physical properties within the entire low temperature regime. The 
Parisi ultrametric structure is shown to determine the dynamic behaviour of fermionic correlations 
(Green's functions) for large times and for the corresponding low energy excitation spectra, which 
is predicted to affect transport properties in metallic (and superconducting) spin glasses. Thus we 
reveal the existence and the detailed form of a number of quantum-dynamical fingerprints of the 
Parisi scheme. These effects, being strongest as T ^ 0, are contrasted with the replica-symmetric 
nature of the critical field theory of quantum spin glass transitions at T = 0, which display only small 
corrections at low T from replica permutation symmetry breaking (RPSB). RPSB-effects moreover 
appear to influence the loci of the ground state transitions at 0(r") and hence the phase diagrams. 
From explicit solutions for arbitrary temperatures we also find a new representation of the zero 
temperature Green's function. This leads to a map of the fermionic (insulating) spin glass solution 
to the local limit solution of a Hubbard model with a random repulsive interaction. This map exists 
at any number of replica symmetry breaking steps K. We obtain the distribution of the Hubbard 
interaction fluctuation and its dependence on the order of RPSB. A generalized mapping between 
metallic spin glass and random U Hubbard model is conjectured. We also suggest that the new 
representation of the Green's function at T = can be used for generalizations to superconductors 
with spin glass phases. Further generalizations of the fermionic Ising spin glass to models with 
additional spin and charge quantum-dynamics occurring in metallic spin glasses, or due to Coulomb 
effects including crossover from 4-state per site to effectively 3-state per site models in the limit of 
infinite repulsive Hubbard coupling are briefly considered. We compare our spin glass results with 
recent d = oo (clean) Hubbard model analyses, paying particular attention to the common role of 
the corresponding Onsager reaction fields. 
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I. INTRODUCTION 



In this paper we emphasize and explore the implementation of Parisi replica permutation symmetry breakingETu in the 
quantum theory of fermionic systems with candom interactions. We present details of the derivation of a many body 
quantum theory with ultrametric structured and we discuss the results in context of the tricritical phase diagram, 
which is given in the subsequent paper II. 

Calculations based on the 4-state per site extension of the standard SK-model provide the framework for detailed 
information that helps to analyze a variety of more complicated fermionic spin- and charge-glass quantum models. 
Spin glasses have been described for a long time by models that isolated glassy magnetic phenomena from other 
physical phenomena. The source of the frustrated magnetic interaction often hidden in the RKKY-interaction of 
metallic systems, this modelling was meant to apply to cases that allow to study glassy magnetic order and electronic 
transport in such a semi-independent way. A first microscopic diagrammatic approach to describe, at the same 
footing the formation of spins and the calculation of their frustrated interaction was designed by HertzQ for a system 
of randomly placed strong Hubbard U centers. A way ofjphenomenological treatment of spin glass order and transport 
theory can be found in the book of Fischer and HertzEL A recent microscopic study of the formation of magnetic 
moments in models of Hertz' type was made by SachdevQ. 

Quantum spin glass models possess a genuine richness in their low temperature phase diagrams that is due to the 
entanglement of a many body interaction and the randomness contained in the couplings. The mixing of spin and 
charge fluctuations leads to new experimentally observable phenomena and to fundamentally new structures of the 
many body theory. 

These systems link phase transitions of remarkably different character: zero temperature so-called quantum phase 
transitions displaying quantum dynamical effects, but also interesting types of thermal tricritical transitionsEJ, which 
mix spin- and charge-fluctuations and also an interference of spin glass features in charge response and transport 
properties. The symmetry classification of these QPT's is totally different from that of thermal phase transitions and 
also does not appear to resemble, for example, the T — Q QPT universality classes of Anderson localization. 
We elaborate details of the manifestation of Parisi symmetry in a T = quantum theoryQ of the fermionic Ising spin 
glass. First extensions to models that possess, beyond the ubiquituous fermion quantum dynamics, also spin- and 
charge-quantum dynamics are also considered. We derive the long-time behaviour of the fermion propagator, which 
proves to be one of the ideal quantities that represent quantum-dynamically the Parisi repiiGa permutation symmetry 
breaking in theories at T = and at low temperatures. The order parameter function (7(a::)lin3 is static and displays its 
nontrivial features only on an interval of order T. It needs the linear susceptibility x — Jq where the spin 

autocorrelation function q remains static for any multi-state Ising spin glass, to incorporate an 0(l)-effect within a 
correlation function. The fermion correlations describe the fermionic spin glasses in a more detailed way than spin- 
and charge-correlations. The latter can however all be obtained from the former ones. Even for spin-static models 
the fermionic correlations, beginning with the fermion propagator, yield a detailed quantttnfcdynamical description of 
RPSB. In addition, this is reminiscent of the classical-dynamical theories of spin glasseaiirEj with imposed Glauber- 
dynamics and despite the fact that the latter was invented to avoid the replica trick thus producing a dynamical image 
of the classical spin glass solution, which indeed represented the Parisi solution. It is this analogy which we emphasize 
by elaborating in this paper the quantum-dynamical image of the Parisi solution for the fermionic Ising spin glass, 
which is spin- and charge-static but has fermion dynamics exclusively determined from the Hamiltonian. Excitation 
spectra encountered in one- and two-fermion correlators are seen to be determined by RPSB in a quantitative and 
quantitative way. 

While most of our results are given for the 4-state fermionic Ising spin glass. The influence of RPSB on the T = 0- 
quantum theory emerges in a general and apparently exists in a model-independent way. This clearly calls for renewed 
efforts to find RPSB in the other domains of disordered electronic systems such as the theory of Anderson-, Anderson- 
Mott-, and similar localization transitions. This is beyond the scope of this paper, but the role of Parisi symmetry 
in the (insulating) fermionic Ising spin glass together with the Ward identity for charge conservation, which we work 
out here, indicates the possible relevance of this symmetry. 

II. MODELS OF RANDOMLY INTERACTING SYSTEMS 

The most fundamental and intensively studied model of random interaction is the Ising model with fully frustrated 
spin interaction. As any spin model it can be described in a grand canonical description with appropriate imaginary 
chemical potential. In this way the standard SK Ising spin glass model is viewed as a fermionic spin glass, given on 
special points of the imaginary axis of the complex /x-plane. On continuing to real chemical potentials one converts 
the Sherrington Kerkpatrick-model into a physical fermionic Ising spin glass. As an interacting disordered fermion 
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model the Hamiltonian appears to be simple. On one hand it only takes into account a random magnetic interaction, 
which, after partial particle-hole transformation is equivalent to a random charge-charge interaction and a magnetic 
field. It allows however for a spin-charge coupling due to redistribution among magnetic and nonmagnetic states, 
which depends strongly on the particle pressure and hence on the effective spin density. 

From the point of view of an interacting disordered fermion system we demonstrate the way in which the known 
instability against replica symmetry breaking leads to a quantum field theory with Parisi symmetry. A large part of 
this paper is devoted to explain this by detailed results for the infinite-range fermionic Ising spin glass. Thus the 
theory presented here starts out from the Hamiltonian, which superficially coincides with the classical SK-model, but 
its spins are defined as the fermionic objects 

(T = n| - nj., ncr = alaa with {a^Q, a^^} = 0, {aj^, 0^73} = SijSap- (1) 

This introduces a new class of observables and correlations which involve odd numbers of fermion operators at partic- 
ular instants of time. Some correlations of groups of even numbers of operators at a given time may become indirectly 
affected, but this type of correlations feels only quantum statistics, which controls the occupation of magnetic states 
I t, >, |0, |> in correspondence with the nonmagnetic ones | t, i>, |0, > and remains static until a noncommuting 
term such as a transverse field or a fermion hopping term etcetera is added to the Ising Hamiltonian. As for the 
fermionic Ising spin glass, the grand canonical description with a chemical potential /i and given by 

JC = H - = - ^ Jtj<Ji<yj - ^ hiai - M X! '^^■'^ 

appears standard, when finally the condition of independent and gaussian distributed random bonds Jij is imposed. 
The physics of this model is however far more complicated firstly thaii^^he 2-state parent model on spin space 
and secondly than related classical spin 1 models, like the BEG-modelta, or the diluted 2-state per site model. 
Complications that arise in understanding the phase diagram of the are deferred to the subsequent paper II. 
According to the statement above, the SK-model and this model prove to be very different even at half-filling. 
Expressing the partition function with the help of Grassman integrals the disorder averaging of the free energy is 
performed by means of the replica trick — [\ogZ]av = lim„^o ^(1 ^ [n^^^Z^Jau). 

The decoupling of eight-fermion correlations generated by the gaussian disorder average is rather standard and briefly 
described in the Appendix. 

The fermionic standard model can be easily generalized or transformed to match more complicated and realistic 
systems. 

While this paper is mainly devoted to the quantum-dynamical representation of Parisi symmetry in the field theory 
of the fermionic rl^i*ig spin glass, earlier studies of spin correlations in more complicated models like metallic spin 
glasses also existlla^ta, and the analysis of fermionic correlations subject to RPSB promises many new results in the 
light of the present results for the insulating cases. 
Writing a generalized grand canonical hamiltonian 

^ = E - E "^S-^^" + ii^lv + SU] E "^T'^^i + E - J)n^.n,,' - E - M E (3) 

ij ij i i,j;a,a' i,a i,(T 

one comprises several subclasses such as itinerant magnetism (with emphasis on itinerant spin glasses), the Hubbard 
model including effects of disorder in U for example. Thus the general model takes care of Coulombic effects, both 
due the Hubbard interaction with and without random fluctuations of the coupling U, and due to the long-range bare 
Coulomb repulsion V(i-j), and secondly allows for transport through either random or nonrandom fermion hopping. 
While the Hubbard interaction for example controls the crossover from 4-state per site to an 3~state per site space 
at [/ = 00 and adds a second magnetic interaction (without frustration if not due to a lattice) , the lopg-range part of 
the Coulomb interaction is also important due to its power to affect or even create gaps for exampleO. Generalization 
due to longitudinal or transverse flelds and due to Heisenberg type of interactions are rather standard and not of 
main interest. Both the metallic Ising spin glass (part 1 and 2) and the Hubbard model are known to be extremely 
hard problems even in the mean field limits and due to 0(1) role of their corresponding Onsager reaction fields. On 
the other hand it is just this common point which, beyond the physical relevance of Coulomb effects in fermionic spin 
glasses—calls for a joint treatment. It is the d = 00-method o£.±he Hubbard model (called dynamic mean field theory 
DMFTlla or local impurity selfconsistent approximation LISAlij by their inventors), which we want to include to some 
extent and compare with the spin glass results. One of the interesting questions concerns the type of competition 
that arises in the limit of vanishing U, which implies an equal weight for local pairing states and for magnetic states. 
The locality of the responsible interaction renders this problem different from frustration in spin glass problems. 
A third group of fermionic spin glass models is given by the Kondo-coupling of itinerant carriers to localized and 
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randomly frozen spins. Since these spins can be represented by pseudofermions coupled to a heat bath by imaginary 
chemical potential, this model falls into the class of interacting systems wi th . tw o fermion species. These s-d models- 
which apply as well to CuMn as to the doped II- VI semiconductor CdMnTeE30, and to heavy fermion systemsE3'EaE3 
are written in tight binding form by 

ij i ij iX ia 

with the localized spins s given in terms of fermion operators c? by = d\-^(T"y di\i where the simplest imaginary 
chemical potential —i^T given here models the Hilbert s pa.ce constraint for spin quantum number i. Higher spin 
quantum numbers can be treated with a little more effortCj'Lj. In the present paper we assume that the localized 
spins are integrated out and end up in a single species model of a type described above. We do not claim that all 
physical phenomena of Kondo systems can be described in this way, but the extraordinary difficulties encountered 
in theory of fermionic spin glasses requires at present the limitation to the simplest models. On the other hand the 
results given here for the fermionic Ising spin glass show that a quantum Parisi phase is to be expected in general. 



III. BASIC ELEMENTS AND SYMMETRIES OF THE MANY BODY THEORY 



Conventional many body theory usually discards effects of statistical fluctuations of interactions such as exchange- 
or Hubbard-interactions. This type of randomness is however not a rare event; a typical exampJfi in the case of 
substitutional disorder is given by the so-called randomly placed U-centers, a picture used by HertaJ in an early and 
pioneering diagrammatic approach to itinerant spin glasses. 

Random interactions can lead to unique and decisive physical properties in all sorts of fermionic quantum spin glasses. 
The entire class of these systems cannot be exhaustively described on the spiin level as the example of effects of complex 
magnetic order on low energy excitations and transport behaviour showedB. 

In this chapter we provide details on technicalities and on physical results obtained by the use of the elements of this 
new type of many body theory. 

One|-p£| the major novelties stems from ther-jaecessity to include broken replica permutation symmetry of Parisi- 
typetrH in the fermionic many body theoryO^Ej. We shall demonstrate that the long-time behaviour of averaged 
fermion propagators is determined by the highly nontrivial Parisi scheme. This bears an interesting relation with 
the time-dependent formalism intended to avoid the replica trick and applied to classical spin glasses with additional 
Glauber dynamics. Our present theory however maintains the replica trick together with quantum dynamics, which 
is enforced by the Hamiltonian rather than being added due to the effect of external degrees of freedom which provide 
a heat bath. 



A. Green's functions and generating functionals for fermionic spin glasses 

The fermionic Green's function, denoted by ^i,cr(e;), and whose standard definition in terms of fermion operators on 
a lattice reads 

G^Ar,r') = -\< Tr{aUr)a],\r')) >1 (5) 

L J ail 

is studied in detail in this paper. It involves apart from the usual quantum statistical average < ... > an additional 
disorder average which uses probability distributions P{Jij) for spin glasses or P{U) for a random U Hubbard inter- 
action for example. In order t be able to study a Parisi type quantum solution, we choose the replica formalism to 
cope with these disorder averages. This formalism builts into fermion field theory the technicalities necessary to deal 
with the typical problems that arise from randomness present in many body interactions. Those are first the presence 
of at least 8-fermion correlations in the replicated action and secondly the broken replica symmetry in connection 
with glassy order. 

The method of time slicing allows to start out from a representation of the fermion propagator in terms of derivatives 
of a generating functional with respect to anticommuting, generating (replicated) fields 77, ^ as 

S3 f 
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where A denotes the effective Grassmannian action. This action is exactly of eigth order in the Grassmann fields, 
since we have chosen a complete gaussian distribution of statistical fluctuations in two-body interactions, 
a) The fermionic Ising spin glass We start with the fermionic Ising spin glass ISGf as our first example, thus avoiding 
complications from quantum-spin dynamics present for example in the fermionic Heisenberg spin glasses or in metallic 
spin glasses of any kind. Quantum spin dynamics, originating from (fermion) spin operators that do not commute 
with the Hamiltonian, are interesting in itself but avoided here to clarify the role of fermion quantum dynamics. The 
latter is shown to be determined on the long time scale by the ultrametric structure of the Parisi solutionO so far only 
known to be present in (static) spin overlap order parameters. We believe that this manifestation of ultrametricity 
in the T = quantum field theory of frustrated interacting systems could also indicate a route to the (to our best 
knowledge) open question of Parisi RPSB in the theory of fermion localization. This refers to the standard part of the 
theory, whereas the question of Anderson localization caused by disorder of frustrated interactions is new and arises 
at a different level; a treatment independent of the precise knowledge of the phase diagram and from the magnetic 
transition would be hazardous. It is clear that this expected phenomenon may be different from Anderson localization 
of noninteracting systems as well as from Anderson-Mott transitions. 

The model of gaussian distributed interaction constants allows to decouple the effective cigth order fermion correlations 
by Q-fields, whose average is static for models without quantum spin dynamics. Fluctuation fields 

5Qf{tX)^Ql\t,t')~<Q-'^> (7) 

are totally irrelevant in this case. 

Taking care of an additional Coulomb interaction (including the local Hubbard part) this situation remains unchanged. 
This will be included below following the effects of a nonrandom Hubbard interaction that plays an interesting role 
even in the absence of a fermion hopping hamiltonian. The cooperation with the nonlocal spin interaction and the 
comparison with recent large d theories of the Hubbard model will then be of particular interest. 



B. Replicated Grassmannian action of the models with gaussian distributed interaction constants 

All models listed in section 2 can be represented by the use of anticommuting fields. The partition function reads 

Z = Jv^^iT)e-^, (8) 

where ^'(t) denotes a master Grassmann fields with imaginary time r generated by time-slicing; this field ^E* contains 
all Grassmann variables that may carry different components in order to take care of different fermion species. Its 
conjugate ^I^ is related to by charge conjugation. 
The action can thus be written as 

I ^ in 

^ = y dT[%lj^^[dr - Hj/ftllp^r + X-yridr " -J^U-yr 

+ H{'lp^^,1p^r,X.yT,X7T)] (9) 

where all "0 and x variables anticommute with each other. The nonhermitean part, which serves to reduce the 4-state 
per site Fock space to a spin ^ space, can be absorbed by U(l) transformation into the fields 

X^{t) = exp{i^^)X',{T). (10) 

In this way the periodicity of the new anticommuting field x becomes 

x{t = /3) = ~ix{T = 0) (11) 

a periodicity that induces semionic imaginary frequencies of the Fourier transformed fields xi^i) with e; = (21 + 
^)TrkBT/h. 

The frustrated magnetic interaction of the fermionic Ising spin glass contributes 

7tT^ Jo Jo 
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to the effective interaction, provided the distribution of the J^'s is gaussian with 2nd moment M{ri — Vj). The spin 
fields are given in terms of the anticommuting fields by 

<M = E ^"aWACaW (13) 

A=±l 



IV. RESULTS OBTAINED FROM REPLICA SYMMETRIC APPROXIMATIONS 

A. Generating functional for the fermionic Ising spin glass 

The derivation of the Green's function is first given in replica symmetric formalism in order to provide the simplest 
possible and Iransparent introduction. The replica-symmetric results are stable solutions only above the Almeida- 
Thouless linfj'Lj and unstable below. The sole presence of random fluctuations in a many body interaction already 
complicates the form of fermion propagators and hence deserves attention before the ordered phase and Parisi RPSB 
in addition will be described. 

We will hence ignore until the next section all sorts of replica symmetry breaking, ie the mostimportant Parisi RPSB, 
also vector replica symmetry breaking (RVSB), recently brought up by Dotsenko and MezardEZi in the different context 
of classical random field systems, which would mean that different copies of the averaged one fermion propagator have 
to be studied. Another type of RPSB could as well arise in form of a nonzero fermion propagator between between 
different replicas. 

The replica symmetric saddle point solution of the Q-field 



Q"'' = g for a 7^ 6, and Q"" = q 
allows to find the Green's function in the form 

5 5 



n/ H' 



,*a('7°,'7°)| 



(14) 



(15) 



where ?/° and z denote replica-local and replica-global spin-decoupling fields and the gaussian integrations over these 
fields are written in the abbreviated notation = -^=e^^^^ . The global field z decouples (Xia VV'^'V'ct)^ ^'^'^ 
^^-integration takes care of VRSB quite naturally within the grassmannian formalism. 
The functional can now be expressed in terms of the Green's function operator g by 



exp 



exp 



exp 



a.l.a 



(16) 



with 



iei +1-1 + a{h + J^z + J^q - qy°-) 



(17) 



For brevity we shall use below the effective field H{z, y°-) = h + Jy^ z + J^/q — q y"". Using the above equations the 
Green's function assumes the form 



/f n„, /!, 9Aei\z,yne^p(.Ea' ^a'{z,y''' \o,o)) 
I^Ujlexpij:a''^a'{z,y'^'\0,0)) ' 



(18) 



which turns into 
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Jy 9a{£i\z,y)exp{<^{z,y\0,0)) 
expmz,y\0,0)) 



(19) 



in the replica limit. 

The rs-functional $(2, y") = $(z, y°|0, 0) is regularized at if = and = in order to retain any physically relevant 
dependence. This results in 



(20) 



The frequency sum evaluated incJ yielded e**^^'^' = 2 
follows from 



cosh((3ii) + cosh{f3H{z, y)) 



which, together with the result 



G rG 



J Z J V 



shows that the replica-symmetric fermion Green's function 
cosh{j3^) + cosh{l3H{z, y)) 



y coshiPfi) + cosh{(3H{z, 0))e5'3'-^'(9-9) ' 



ri^i\z,y) 



(21) 



1. Fermion Green's Function in the disordered phase (T > Tf) 



The fermionic Ising spin glass does not allow for spin dynamics unless Glauber dynamics introduces it by coupling 
to a heat bath. This would be analogous to the SK-model. The JS'G/'s Green's function however is always dynamic 
and the way in which dynamical behaviour shows up in the spin glass is however nontrivial. The simplest task is 
to determine first the Greens function above the freezing temperature and to compare it for example with that of a 
random Hubbard model in the local limit. The fermionic Ising spin glass allows for an exact evaluation of the bare 



fermion propagator C/j, 



< 



>]av in the disordered phase; here H denotes the usual effective field and the 



double average refers to the replica-local and the Parisi block decoupling fields. The result can be written in the form 



Qij.cr (^n) 



(22) 



(2 - A2)e^(^-^^)'g~c/i((l - X^)P^ia) 



ifia + iXpj^q (- 
-))e 



2qJ 



^.1 g 



where fi^ = 11 + aH includes a magnetic field H and e„ = (2n+ VjirksT /Ti. By analytical continuation to real energies 
(the error function with imaginary argument may then be expressed by the conflueB|t hypergeometric function with 
real argument) one easily extracts the disorder averaged electronic density of statesE^ 



pAe)^--ImG'^{E) = 
and the real part of the Green's function as 



ch{(3^j.) +ch{P{e + fi)) 
ch{(3fi) + ch{(3H)exp{^(3^.Pq) 



(23) 



TZeGaiE) 



ch{f}^i„) + ei0'J''^ '2' 2' 2J2g 



expap^J^q- 



)E + (3.Pq ^,13 {E + p.Pqf 



ch{Pn„) + exp{\(3'^,Pq) 
exp{\f?JH-^^^iff-) 



2Pq 

E - liPq 



chiPfi^) + exp{^(3'^,Pq) 2Pq 
where IFl denotes the confluent hypergeometric function. 



iFii 



2' 2' 2Pq 
1 3 {E-(iPqf 



2' 2' 



2.Pq 



(24) 
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S. Spin- and charge-response function derived from the spin glass fermion propagator 



As a basis for calculations in the spin glass phase, and in particular under respect of replica symmetry broken fermion 
propagators, we demonstrate the easiest example of spin- and charge-susceptibilities above the freezing temperature. 
Both correlations given by 

Xc EE T'iY^ 9l{ei)Uiy) - [{Tj2<T'^9a{ei) fU(y) (25) 

cr,€ 

= T'J2^liei)-{Tj2a<g{e)f (26) 

cr,e ere 

are static and satisfy the relation Xa- + Xp — '^'^ = 2[h]. 



3. Fermion Green's function in the spin glass ordered phase (T < Tf) (replica-symmetric approximation) 

Despite the fact that the replica symmetric result is unstable below the fermionic Almeida Thouless line, hence in any 
case below the freezing temperature, the fermion propagator is complicated to an extent that requires insight from 
its simplest realization. This is the extension of Eq. ( p2| ) by including the effect of one replica symmetric spin glass 
order parameter q. By evaluating Eq.(|^) we obtain on the imaginary frequencies 

t / , _ 1 ,2 i 



4 J J-oo ch{f3fi) + Exp{\p^J^{q - q))ch{(3J^z) 

(.;-i(^.+jv?z))2 ei - i(n+ J,/qz)^^ 

J^/2{q - q) 

The analytic continuation of this expression provides us with the replica symmetric density of states, showing a 
magnetic hardgap as discussed below, and the real part of G^"^ (E) given by 

Re\G^''\E)] = ^ / dzExpi^-z^)^— t^-t^t r-— r (28) 

^ V2^.J{q-q) J^oo 2 ' ch{(3y)+Exp{}^l3^J^{q-q))ch{(3J^z) ^ ' 



Exp 



\m r ^^R^l-^(~ \ {E + J^z + XpJ{q-q)) 
Xf3Jy/qz + -P J (q-q) 



2" 2J^{q-q) 



\=±l 

13 {E + J^z + XpJ{q-q)f 
'^4' 2' 2j2(g-9) >> 



The asymptotic behaviour either for large energies or for T — > can be checked by the use of 1 — Er f (x) ^ --^2££LjlJ. 
and the functional relation Erf{x) — -^1^1(5, |, — x^)!!. 

The result Eq.(|2|) is displayed in Fig.(0). Comparison with the extrapolation of the (paramagnetic) q = 0-result 
below Tc, the role of the (replica-symmetric) susceptibility as a measure of the spread of just two wells forming below 
the freezing temperature is shown. This spread is seen to be related to the size of the magnetic hardgap in this 
approximation. We shall come back to the essential changes of the fermion propagator in the small energy- and in 
the long-time regime below in the section on replica permutation symmetry breaking and Parisi symmetry. 
The fermion density of states in replica symmetric approximation already reveals the strong coupling between spin 
glass order and fermionic properties. The strongest effect emerges in the zero temperature limit. After analytic 
continuation je/ — *■ e + «0 we find 

coshiPn) + cosh{f3E) f°° e 



P<y(e) = — / dz !„■>,■>,- — 7 (29) 

2TTjy/q^ cosh{(3^i)+cosh{(3H{z,0))ei'^'J"('i-i^ 

The T ^ limit is not simple. We first consider the regime 
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\p\<lf3J^{q~q)- (30) 

One should note that q — q — 0{T), whence this condition requires the Fermi energy to satisfy 



This value was obtained from the rs-selfconsistent order parameter equations above. One may prefer to identify 
l3J{q — q) with the linear susceptibility x of the ferniionic Ising spin glass ISG / . 



e-iP^-^''^^-'^^cosh(f]n) +cosh((3E)) , 
P^^'^ = 2nJVq—q I 



cosh{l3H{z,0)) 



(32) 



The T ^ limit is most easily obtained by integrating over the rescaled variable u = f]H{z,0). This isolates the 
leading contribution as 

cosh(Bu)+cosh(BE) _ i (fl^ /a^^ , >, 

p.(e) = rr^—^ " «'+^+^^V(r,i?), (33) 

where the integral 

nT.m=r,u ,34) 

"'-oo cosh{u)e ■''■'(9-1) 

converges for energies satisfying the condition 

E<l3J^{q-q). (35) 

Since the maximal \E\ within this regime coincides with the limit of a vanishing prefactor and hence a vanishing 
density of states in the T — > limit, this energy is to be identified with the gap edge. We call the gap energy Eg{h). 
The replica-symmetric solution for energies inside the T = 0-gap (and Fermi energies within half of this gap) is hence 

EgW := Je-^^ (36) 

and the final result for the leading low temperature DoS-contribution, which vanishes exponentially within the gap 
\E\<Eg[H), reads 

Two striking aspects emerge: 

i) there is neither a spin dependence of the gapwidth nor of the leading and exponentially decaying low temperature 
density of states despite the magnetic field dependences; 

ii) the given result is valid only for < ^Eg{H), while there is no limitation on = |e + Thus the conclusion 
on the gapwidth is limited to a Fermi energy varying over half the gapwidth. 

iii) one would like to see how the Fermi energy can be moved through a gap edge in order to obtain a non half- 
filled ground state. From the analysis of the free energy we know however that phase separation occurs beyond 
|/xo| = ^Eg{H), at least in the replica-symmetric approximation. Hence we believe that instanton solutions may have 
to be used, replacing the homogeneous saddle point solutions above. We conjecture that an instability, identified as a 
second negative eigenvalue of the Hessian matrix - the first being just the known Parisi-RSB instability -, indicates 
just this breakdown of the spatially constant solution in the regime of first order transition from a half-filled to a 
completely filled (or empty respectively) system. 

One should nevertheless realize that only ^ = corresponds to half-filling at any temperature; moreover it is interesting 
to compare the domain of continuous thermal spin glass transitions which extends to the tricritical points at either 
HTc3)\ ~ -96125 at T^s = |J with the \^iiT = 0)| = « 0.398942J. 

Real and imaginary parts of the Green's function are plotted for various cases in Figs.(|l|,^,||). 
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Re G 




FIG. 2. Cross sections through preceding figure (rotated) of Re[G] at temperatures T = oo (curve a), T — Tf (b), T = .5 J 
(c), and T = .2 J (d) 
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FIG. 3. The imaginary part of the replica-symmetric (~ DoS) as a function of energy S = e (/i = 0) and (5 J = ^. 



p(E,H = 0) 




-4-2 2 4 

E 

FIG. 4. Creation of the spin glass related gap as the temperature falls below T/ (cross sections of preceding figure); curve a: 
T = oo, b: T fa T/, c: T = J/2, A: T = J/3, f: T = J/6, T = .IJ, g: T = .05J , and h: T = .025J. 

4- The magnetic hardgap of the DoS at T — 

The density of states p{T — 0, \E\ > Egih)) is nonzero; it can be evaluated by the saddle point method. We start out 
from Eq.(|29|). The very same condition < ^Eg{h) applies again. We shall come back to this below. 
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While q, x, and can be obtained by extremizing the free energy, ie without exphcit knowledge of the underlying 
density of states (DoS), this fermionic DoS is of course related to the magnetic properties of the system and its 
analysis requires the solution of several selfconsistency equations for spin- and for charge correlations. We obtained 
these solutions for low temperatures and for the entire range of magnetic fields and of single fermion excitation 
energies. The existence of a hardgap of size 2Eg{H) is best appreciated in our T = 0-result for the density of states, 
which is given by 



1 



27rJ 



Q{\E\ - Eg{H))exp 



\E\ 



(38) 



with E :— e + ij. For finite range models we expect exponentially small corrections within the gap and due to Griffiths 
singularities. Eq.(|38|) is valid within the regime < ^Eg, which corresponds to a half filled system at T = with 
filling factor 



z/ = l + (l — g — Tx)tanh{(3ii). 



(39) 



The zero field ratio Eg{0)/Tc increases with the chemical potential from 1.179 at /i = to 1.238 at ii — J/v27r. In 
terms of Eg{H) the following low temperature expansions (for linear susceptibility and order parameter), required for 
the exact evaluation of the DoS formula (E9f) at low T, were obtained as 



J J 



y2 



X 



J2 



2 J3 



(1 



0{ 



J2' 



(40) 
(41) 



These results show that i) the density of states is zero at T = in the finite interval given by |e + /x| < Eg{H) and 
ii) the gapwidth shrinks as the magnetic field is increased. This will be the source of the negative magnetoresistance 
in the extended Ising spin glass model with charge transport as discussed below. The physical consequences of these 
spin glass related properties are hence in agreement with the experimental observations of crossover behavior in the 
low T resistivities mentioned above. For single fermion energies E = e + fj, smaller than the gap energy Eg{H) and 
< ^Eg{H) the density of states decays to zero exponentially as given by 



P^(e) = 



1 



2 Jcos(| 



E 



■{/3Eg{H))-^ [ch{f3fi) + ch{PE)] e 



(42) 



2 Eg{H) ' 



We find that the hardgap persists also in itinerant models with a fermion hopping term added to the nonconducting 
Ising Hamiltonian, until the bandwidth exceeds a critical value. The gap and its related properties of this iSjctepded 
itinerant spin glass model, whose magnetic phase diagram and -transitions have previously been analysedljM, are 
discussed below. 
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E 

FIG. 5. Single particle density of states (DoS) for the fermionic Ising spin glass {ISGf) versus energy and magnetic field in 
units of J. The hardgap centered around zero energy is due to random magnetic correlations and manifests itself in an activated 
behavior of the low temperature hopping conductivity of disordered magnetic materials. The reduction of the gapwidth by an 
external magnetic field explains the observed negative magnetoresistance. 

Coulomb interaction effects are of diverse nature: the long range part will still tend to depress the (remaining) DoS 
near the Fermi level and thus stay particularly relevant when either spin glass order is weak or absent, or (in the case 
of fully developed spin glass order) when the Fermi energy lies close to the gap edges. The Hubbard coupling U leads 
to a shift of the chemical potential {fj. ^ fi — ^-i^), and also to a shift H ^ H' = H + !^m{H) of the applied magnetic 
field. This implies for finite temperatures that the Almeida Thouless line depends on the Hubbard coupling too. 



5. The replica-symmetric approximation of the T = fermion propagator 

The density of states could be used as the spectral weight in the Kramers-Kronig relation for the fermion Green's 
function. One may as well extract the zero temperature limit of Eq.(p7|). This results in the (bare K — 0) fermion 

propagator Oi^'^'' (e) 



dz _ 1 ^2 
=e ^ 



{- 



i{fi + z + xo) e-i{^i-z-xo) 



} 



(43) 



with energies e; = {21 + 1)'!tT — > e lying dense on the imaginary axis and xo = limT^o(3{q — q) representing the 
RS-approximation of the susceptibility (here and in the following we set J = 1 for brevity) . The most remarkable fact 
being that the gaussian distribution of say self-energies z ~\- xo are now cutoff, which is responsible for the hardgap 
described above. Recall that in Eq.(|2^) one starts off from a comple te g aussian average over all z. 
Insertion the zero temperature results for the density of states, Eqs.(29) and (^, in the spectral representation 



de 



iei - e' 



(44) 



confirms the representation of Eq.(^). The final integral in Eq.(|4^) can be evaluated to give the following explicit 
solution 



1 



E 

A=±l 



inXsign{e)e' 



(45) 
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with 



Ax{e) = €-l[^l + Xxa\. (46) 

and Ei(x) denoting the exponential integral functionil. The results are given for zero magnetic field; the fact that the 
chemical potential appears in the form /i ± xo shows the effect of the Onsager reaction field. The oversimplification 
of the replica-symmetric approximation is revealed by the /i-shift, resembling only contributions from up and down 
local fields. 

The time-dependence is obtained from 



dz 



27r 



+ [e-^+"0(^+(z)) + e- 
where we introduced effective chemical potentials 



/i±(z) = ^± (z + Xo) 



(47) 



(48) 



that depend on the spin decoupling field and on the lowest order susceptibility. The latter represents the Onsager 
field contribution. 

We recall that the given solution is valid for ^ < ^xo (apart from its instability against RPSB that is taken into 
account below), which means half-filling at T = 0. Evaluation of the imaginary-time dependent propagator yields 



1 



-(A'-Xo)t-^ 



1 + Erf{ 



T2' 



(r) + le 



-(A'+Xo)t-^ 



1 - Erf{ 



V2' 



9{-t) 



(49) 



Alteration of this time-dedendent behaviour under RPSB will be an important conclusion below. Before this problem 
is attacked, it is interesting to compare the replica symmetric fermionic spin glass propagator with a propagator of 
random U Hubbard models in the local limit. Fermion hopping, which introduces the full complexity of the Hubbard 
model, plays also an important role in promoting the insulating fermionic spin glass to a metallic spin glass model. 
Since all of the complicated physics of both model classes is sandwiched between the (identical and solvable) fermion 
hopping model and the (different and solvable, provided an infinite-range magnetic interaction is considered) pure 
interaction limits, the latter limit offers the most interesting case for comparing the two solutions. 



6. Comparison and mapping with the Green's function of a random U local Hubbard limit 

The wellknown exact SQ luti na. of the local limit played an important role as starting paints of recent many body theories 
of the Hubbard modcHH^cJ. The field-theoretic decoupling given by VoUhardttS can be extended to account for 
disorder fluctuations in the coupling U. We shall compare the spin glass and the random U Hubbard model field 
theoretic techniques in Appendix B. For the purpose of comparing the replica-symmetric fermion propagator of the 
fermionic spin glass with a potential random U Hubbard analog in the local limit, we do not need this apparatus. It 

~ As 
imit. 



is sufficient to consider a gaussian average over U with a properly chosen cutoff to see the similarity with Eq 
well one may compare the Green's function obtained above the freezing temperature with the local Hubbard 
At half-filling the Green's function reads 



GAei)^k . A^/. + ■ \rJ (50) 
2 lei + u/2 lei — (7/2 

Since this is the exact solution for any half-filled U-realization in the local random U model, it is easy to evaluate not 
only gaussian distributions but also Lorentz-, box-, and semicircular-distributions of random U-couplings in order 
to obtain the disorder-averaged one particle Green's function. The average < ^017.^ > can also be derived. Thus these 
results serve also as a interesting test field for the replicated field theory, which usually runs into formal difficulties 
whenever nongaussian distributions are involved. 

For unrestricted gaussian-distributed U (equal weight of attractive and repulsive interaction) one obtains 

[gA^i)]p^(U)=^e^'/''"il~Erfi^))sgnieO (51) 
and the density of states is a simple gaussian 
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\p]p.iu){^-\j^^-''^"'" (52) 

In order to generate a solution that matches the fermionic spin glass propagator, Eg. (^) , one needs to cut off the 
gaussian distribution. i_. 

For /i = it is just the Onsager reaction fieldEj of the spin glass, which requires to introduce a finite average < U > 
by 

+SU, (53) 
where the fluctuations obey a gaussian distribution of nonnegative (5C/-fluctuations 

P„^((5C/) = -^exp{-l{6Uf)e{6U) (54) 



V ^TT ^ 

which maps the fermion propagator of this random U Hubbard interaction onto the replica-symmetric approximation 
of the fermionic Ising spin glass propagator at T = and /i = 0, provided (again given in dimensionless quantities) 

<U>^ 2x0 (55) 

This mapping involves only repulsive but random interactions. Small random deviations from the average positive U 
are strongly weighted, which renders the model perhaps a realistic one. The relation (|5^) also provides an illustration 
in terms of the Hubbard interaction for the gap-width 2xo of the fermionic spin glass. 

It is interesting to imagine an additional fermion hopping term, thus a Hubbard model with the given U-fluct nations: 
realizations of large U can be described by a t-J model and the entire regime down to perturbative small U is contained 
with increasing weight. 

The random U local Hubbard limit, introduced here for purely formal reasons, is a zero-dimensional problem and 
one cannot expect to find a spin glass transition like the one described for the Ising spin glass model with nonlocal or 
even infinite-range interaction. Allowance for fermion transport however reopens this possibility and requires further 
analysis, also in relation with other metallic spin glasses. Moreover, the problem of RPSB can be raised for the local 
limit of the random U Hubbard model too, which means that the exact solution in the random case may not simple. 
We defer this analysis to another publication. 

In contrast to the described situation, which contains a density of states gap identical to the fermionic spin glass by 
construction, a gaussian distribution without cutoff removes the gap and poses hence a problem to the existence of 
a Mott transition of the Hubbard model. A striking difference with respect to the interaction gap of the frustrated 
nonlocal spin interaction is to be noticed. In this respect the statistical fiuctuations of the Hubbard interaction 
compete with both the nonrandom part of U and with the spin coupling Jy . We first need to evaluate exactly the 
partition function and the generating functional of the random local limit, which passes through a time-dependent 
Grassmann field theory. 

In contrast to fermionic spin glass models it is easy to consider distributions other than gaussian in the local Hubbard 

limit. The latter may thus become helpful for the spin glass problem. 

For example, the box distribution Pb{U) yields the averaged Green's function 



U2-U1 

which turns into the retarded function 



U2 Ui ' 

arctan arctan — 

2ei 2ei 



2e - C/2 2e - Ui 



2e - C/2 2e + Ui 



(56) 



(57) 



Depending on the parameters of the box distribution the logarithmic cuts may either overlap or remain separated. In 
the latter case the Mott transition induced by sufficiently strong fermion hopping will occur; the case of a local limit 
with arbitrarily small gap should help to study the Mott transition also arbitrarily close to the local limit. 

7. Some physical quantities related to the gapped fermion density of states for T < Tf 

Below the freezing temperature 0((7^)-corrections occur in < p >. Using the exact low temperature solutions of 
section 2.2 one can obtain some information about the density of states deep in the ordered phase. From the relation 
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[i^]av = 2 — g (valid in this form only at T = 0) and q{T = 0,0 < ^ < -^=) = follows that the system stays half filled 
independent of the chemical potential. This quantity can also be obtained from the disorder averaged one particle 
density of states by 

/C30 
def{e)[pU{e,y) (58) 

<y=\,i 

and hence display consequences of the interaction induced gap around zero energy. j— . 
The situation may be comparable to the antiferromagnetic gap in the clean Hubbard modelEj. 
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V. QUANTUM FIELD THEORY WITH BROKEN REPLICA PERMUTATION SYMMETRY 



This chapter describes the central results of the present paper. The Parisi symmetry, described inO in terms of 
symmetric groups and more commonly known for its ultrametric structure, is now shown to characterize a quantum 
field theory. The problems which arise in the QFT due to the effect of RPSB are first resolved in a one-step breaking 
scheme. Then general relations for arbitrary steps of RPSB are found in addition, which allow conclusions on the low 
temperature behaviour. 

The fermion Green's function can be derived, like any other correlation function and at any order K of Parisi-RPSB, 
from the generating functional S. The one-fermion Green's function can be obtained according to 

^M>(^0 = -4l7i'"S„(M,{r7}) (59) 

employing derivatives with respect to the generating fields rj and 77. This functional involves gaussian averages over 
(K+1) fluctuation fields z-y. Since details are given in Appendix 1, we simply introduce its structure by the following 
symbolic shorthand notation 

S„(r;,,y) = e-^''^^^^'^«-- 

where 

5.(p, ei\{z-,}) - {iei + + <7H[{z^-})r^ (61) 

denotes the unaveraged fermion Greens function of the localized fermion exposed to the statistically fluctuating ef- 
fective fleld 



Jl ^ j V^Exp [^'g"^^' + J/^- - ^f] (60) 



H{{z'^}) = H + Y,^q,~ g^+i (62) 
7 

which depends on all fluctuation flelds z^. The fermion-propagator of the insulating models is of course local in real 
space. Spin (decoupling)-fields z^, carrying a Parisi block index 7 (a^ runs over all replicas which belong to a particular 
Parisi block), explore the random magnetic order. In the metallic case, the bare Green's function go = g{H = 0) 
and hence g itself become nonlocal, a complication that need not be studied now in order to understand the Parisi 
symmetry of the QFT. „ 

The Parisi matrix Q parisi has the wellknown forma apart from the nonvanishing diagonal elements q (which we identify 
with go); their presence is required by the fact that (ct^)^ = (n| — fii)^ ^ 1. The structure of the Parisi-matrix is of 
course responsible for the rather complicated form of the Lagrangian (see Appendix 1); despite this complication the 
fermion fields can be eliminated in the standard way, which leads to the selfconsistent equations given below. 



A. One— step replica permutation symmetry breaking 

It is instructive to start with a one step replica symmetry breaking (denoted by K = 1), which is also considered 
to be the correct result in the regime of discontinuous thermal transitions. Here we wish to develop the full Parisi 
solution from this first nontrivial starting point. Denoting by qi and by q2 the new offdiagonal block elements, the 
decoupling field z splits up into zi and Z2. The Green's function assumes the suggestive form 

g„{ei) = [[[{iei + + aH{z^,Z2,yT^)]yAz^U (63) 

which includes three successive averages to be specified below. The instability of the replica-symmetric solutions be- 
comes strongest as T ^ 0, but only recently this has found its explicit representation in a zero temperature quantum 
field theory including the low temperature effects. While quantum magnetic transitions are so far unaffected by this 
instability to a large extent, the fermionic density of states and hence the fermion Green's function were observed to 
depend strongly on RPSB. This leads to a quantum field theory with Parisi symmetry, ie with ultrametric structure. 
The fermion propagator had acquired a nonsimple form already in the replica-symmetric approximation due to the 
random interaction on one hand and furthermore through the effects of spin glass order. We now wish to show in 
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more detail how this solution is improved at first-step breaking-iii' =1) ^^^id furthermore how a K~invariant relation 
can be derived to prove the one-fermion pseudogap reported irJj. 

We have recently presented the first solution to the question how Parisi replica permutation svmtoetry breaking 
(RPSB) and the related nonconstant part of the Parisi spin glass order parameter function (j(a;)0'ErH are displayed 
in the low temperature many body theory of fcrmionic systems with frustrated Ising-interactions, emphasizing the 
T = 0-limit in particular. Despite the fact that the interval < x < xi oi nonconstant q{x) representing RPSB 
vanishes with temperature T as T — > 0, we find a large 0(r'')-effect to persist in many important physical quantities. 
This includes replica-diagonal fermion Green's function and fermion density of states, where at any step K of RPSB 
the set of different order parameters is seen to determine the quantum-dynamical behaviour of the fermion propagator 
and of vertex functions. These effects are complemervtary to and not ia contradiction with recent replica-symmetric 
descriptions of T = quantum spin glass transitionslij. Parisi-RPSBa is seen to decide the qualitative and quanti- 
tative features of the low energy excitation spectrum. While results are presented for an insulating model, the effect 
appears to be rather model-independent and should hence be felt in transport properties of models with additional 
hopping hamiltonian for example. 

We wish to provide results which evidence the fact that fcrmionic spin glasses also link closely glassy magnetic order 
and transport behaviour; further similarities between Hubbard model ajod-the fcrmionic spin glass have been traced 
back to the particular role of the Onsager-Brout-Thomas reaction fieldallj for all these systems, as can be observed 
by comparing Hubbard-CPA-E3 with fcrmionic TAP-equationstll. 

Spin- and charge-excitation spectra of fcrmionic spin glasses must be evaluated in order to construct a meaningful 
many body theory. This article focusses on the effect of Parisi replica permutation symmetry breaking (RPSB) on the 
single fermion density of states (DoS) , hence on the fcrmionic Green's function, and, by virtue of the Ward identity 
for charge conservation, also on vertex functions, thus on the entire ensemble of quantities that provide the basis of 
many body theories for fcrmionic systems with frustrated interactions. 

It is known that, jijeplica-diagonal quantities like the linear equilibrium susceptibility x feel Parisi symmetry breaking 
even at T = Cuij'EI despite the fact that the nontrivial part of the Parisi function only lives on an interval of width T. 
The susceptibility had been analysed by Parisi for the standard SK-model who found a rapid convergence towards 
the exact result as the number of order parameters increased, this number being equal to if + 1 in the SK- and 
equal to if + 2 in fcrmionic models. While the low temperature regime of the SK-model had not been of particular 
interest from the point of view of phase transition theory, it becomes highly important for fcrmionic spin glasses, since 
the T — 0-theory of excitation spectra plays a crucial role and, for the additional reason that some models exhibit 
quantum phase transitions along the T = 0-axis. Parisi nevertheless analysed the low T regimetl of the SK-model 
finding that K-step RPSB on one hand provided increasingly good approximations but failed to completely remove 
the negative entropy and the instability problem at low enough temperatures unless K oo. 

In this paper the effect of one step RPSB on the density of states is presented in detail, followed by the derivation 
of an analytical relation valid for all K, which allows to determine the type of excitation spectrum present in the 
full Parisi solution for the fcrmionic Ising spin glass. Despite the fact that the regime of deviation from a replica- 
symmetric spin glass order parameter is only of 0{T), we find that it has a large 0(T'' )-effect on the density of 
states, the one-particle-, and many-particle Greens functions at T = 0. This density of states is derived as usual 
from the imaginary-time (disorder-averaged) fermion Green's function [< ■0(T)^(r') >]av, which is one of the decisive 
quantum-dynamical elements of any many body theory of fcrmionic spin glasses. This illustrates that, unlike the 
usual picture of a Parisi solution being just a static order parameter function, the fcrmionic picture must include the 
qualitative extension to dynamical quantities. Those become drastically altered by the nontrivial part of the Parisi 
solution, which is otherwise invisible at T = 0, hence providing a quantum-dynamical image of RPSB. 
Let us consider the Parisi solution of the infinite-range fcrmionic Ising spin glass defined as the SK-model extension 
in a grand canonical ensemble with spin one half operators represented as ct^ = fi| — . 



1 . Fermionic density of states 

It is known since Parisi's work! that an analytical low temperature expansion is hard to obtain even for the standard 
SK-model and its smaller set of selfconsistcnt parameters. First insight is gained by the one-step RPSB {K — 1). 
The standard three parameter set of the SK-model for K = 1, order parameters qi and q2, and m = mi ^ T, is 
enlarged in the fermionic space by g — gi T, where q := [< (T(r)cr(T') >]av represents a spin correlation, which 
remains static unless a fermion hopping mechanism or other noncommuting parts are included in the Hamiltonian. 
For the fermionic Ising spin glass the {K = l)-DoS reads 
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ch{0fi) + ch{(3E) J (9-91) 



^2TT{q-qi) 



dv2e ^"2 



J —00 ^ 



where 



C = cosh{(3H) + C, 



C = cosh{P^i)exp{--0^{q - qi)). 



(64) 



(65) 



This expression reveals the competition between the particle "pressure" exerted by the chemical potential /i and the 
single-valley susceptibility 



(66) 



leading to a crossover at — iii the T 0-limit. The C~term is a fermionic feature which is absent from 
the standar|d]SK-model. It is closely related to the fermion filling; this filling factor behaves discontinuously on the 
T = O-axistj. The zero temperature limit allows to perform the wi-integrations, which results in the exact T = 
formula for the density of states 



^-ia^(ff)(l-g2)- 



using the definitions 



1-92 



7:^1 -q2{H) 



— i ^ ( 

3 2 1 — Qo W — g 



dz- 



d{z) + d{-z) 



-e(Ai 



d{z) 



1 + Er}[ 



a{H){\-q2) 



922, 



^2(1 - 92) 



Ab = 1^1 



d 



a{H) = lim^_^„—m{T,H). 



(67) 

(68) 

(69) 
(70) 



These equations allow to derive the first improved approximation beyond the replica-symmetric solution, given in 
the preceding section, which displays a magnetic hardgap of width 2Eg(H) in the DoS and the fact that half-filling 
at r = extends over the finite interval of chemical potentials given by < ^Eg{H). Beyond this interval phase 
separation occurs together with a discontinuous transition into a full or an empty systemE£l. More details about the 
phase diagram are presented in the subsequent paper II. 

A stable and spatially-homogeneous saddle-point solution could only be found for the half-filled case. Thus the 
following analysis is restricted to this interval of chemical potentials. Its width is determined selfconsistently and seen 
to decrease to zero as K tends to infinity. For < i limT^oiP{q — qi)) '^6 first derive the coupled set of selfconsistent 
equations 



= 1 - 



91 



92 = 







dr 



-f 



The shorthand gaussian integral notation 



{qi-qi) + 



duh 



T 



In / C 



T 
m 



C"HnC 



^2TT{qk - qk+i) 



exp 



2{qk - qk+i) 



(71) 
(72) 

(73) 
(74) 

(75) 



adopts the normalization used by Parisiu. For T ~ 0, we solve the selfconsistency equations up to one final integral, 
finding the T — 0-set of coupled equations 
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q = qi = l, lirriT^o 



Q-Qi 



= 1- 



T 

°° dz 



= X, 12 = 



°° dz _<f2fV!i2! 

e 2 



d{z) - d{-z) 
d{z)+d{-z) 



oo 



27r 



{--ln[-e-^'^\d{z)+d{-z))] 



+ 



{at + ^ z)d{z) + {at-^ z)d{-z) + ^/stj^ g-ia^t-l 9,^7*1 /[d{z) + d{-z)]} 



where t = qi — g'2-All parameters represent the temperature- and magnetic field-dependent solutions. 
The solutions for T = and H = are given by 

52 = 0.476875 

a = limT^o-^m{T) = 1.36104 
w = limT^o^^Y^ ^ .239449. 




g.4T = 0,H) 



iii(T = 0,-H) 



H/J 

FIG. 6. Field dependence of dm/dT (top), of the order parameter 52, and of gapwidth parameter w (bottom) for 
and zero temperature. 
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E/J 

FIG. 7. Zero temperature density of states as a function of energy and magnetic field for 1-step RPSB 

The H-dependent solutions shown in fig. 6 are then used in evaluating eq.(^^ for the density of states. T — 0-results 
are shown in figs. 7 and 8, while the result at finite low temperature of fig. 9 illustrates the presence of plateaus of 
constant slope, each corresponding to Parisi order parameter separations (here: q — Qi and qi — (72)- The number 
of these plateaus of constant slope increases with the order K of Parisi-RPSB. Hence, the time-dependence of the 
Green's function should characteristically depend on the order parameter separations qk — qk-i- 
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FIG. 9. The 1-step RPSB solution for the zero field density of states at very low temperatures (T = .OIJ) (p is symmetric 
w.r.t. E). The T = hardgap of reduced size is already digged out, yet two different slopes ^ corresponding to qo — qi = 0{T) 
(steep portion) and to qi — q2 = 0{1) remain visible 




-4-2 2 4 



FIG. 10. DoS at l^tep RPSB for magnetic fields H/J = .2 (curve a), .3 (b), .5 (c), and .6 (d) 

If we compare with the repUca-symmetric result a reduction of the gapwidth is observed. Analytically one finds a 
gapwidth 

Eg{H) =x = UrriT^oPiq ' 9i) (81) 

which turns into UmT^o{(3{q — qiX))) in terms of the Parisi function q{x) B.t K = oo. Only in absence of RPSB this 
susceptibility coincides with the equilibrium one, denoted by x- In f^ct for 1 step RPSB the fermionic Ising spin glass 
approaches x = /^(s — 1i) + PnT,{qi — (72) — * -95, the same numerical value as the one for the SK-model. 



2. The fermion propagator 



We are interested in the explicit form of the fermion propagator at one-step RPSB for several reasons. To mention a 
few, this function Q~^{ei) is an indispensable element of the quantum field theory with Parisi symmetry, it is directly 
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linked to some observables, it contributes by Re(G) a second quantity which illuminates the crossover to a pseudogap, 
it is also needed to formulate the Ward identity, and it add s another challenge for the mapping on a random U 
Hubbard model in the local limit. We start out from Eq.(124) (as derived in Appendix 1) which yields in the replica 
limit for K = 1 



where ea;p(4>) = 2{cosh{f3H) + cosh{(3^)e 2 
In the regime /i < one obtains 



2 



(82) 



J zi 



A=±l 



1 - Erfi 



e, -i{fi + H„ + Ax) . 
v/2(9-9j ■ 



(83) 



The zero temperature limit (using limT-^oPm — a finite) simplifies this result and can be cast into the form 



J— 00 ^ 



(84) 



where Hq = H{0, zi, 22)- The spectral representation ( ^ ) now employs the f-step RPSB-result for the density of 
states Eq.(||,|7|). 

Before evaluating the zi-integrals we can reconsider the mapping with the local limit of a random U Hubbard model 
and the changes that occur due to RPSB in the spin glass. Now 



<U > 



X 



(85) 
(86) 



The second correspondence shows that only positive [/-fluctuations around a positive mean value occur, which de- 
creases with each step of RPSB like the gapwidth of the fermionic spin glass. However, modelling the latter by the 
random U interaction requires {K + 1) fluctuation-fields Sua at K-th order RPSB and subject to the constraint 
6U = ^"^a > 0. Apart from the overall condition of positive U and SU the fluctuation-fields are allowed to assume 
any negative value. Let us express the speculation that the existence of realizations of arbitrarily negative attractive 
interactions might open the way to a pairing creation due to Parisi symmetry breaking. These pairs might become 
delocalized due to the introduction of a fermion hopping term. We do not pretend to have an answer at the moment, 
but we think it is justified to raise the question whether a superconducting transition might become possible due to the 
Parisi replica symmetry breaking at and beyond the instability line towards RPSB. 

The random U local Hubbard limit was constructed to match the fermionic spin glass propagator at T = 0. It is clear 
that this mapping cannot be achieved selfconsistently for all temperatures, since the zero dimensional local limit does 
not allow for a thermal phase transition. We consider it possibly important that one may achieve a mapping at zero 
temperature between a random U Hubbard model and a metallic Ising spin glass with identical hopping hamiltonians. 
This will involve the question of replica symmetry breaking in the random U Hubbard model. 

The role of frustration and the analog of the Almeida Thouless instability in the random U Hubbard model arc to be 
studied, which requires to study the replicated field theory with incomplete or nongaussian disorder average. This is 
a problem for later studies. 

We also consider interesting the question, whether frustration in clean systems leads to time-dependent behaviour 
that resembles the one generated by replica symmetry breaking. 



3. Quantum-dynamical image of RPSB 

The results for low energy excitations indicate that the time-dependence of fermion correlation functions must be 
affected by replica symmetry breaking. We study now the fermion propagator and in particular the retarded Green's 
fucntion. The Fourier transformation 
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G^W = / ^e-*^*G^(e) (87) 
-oo 27r 



yields 



W^g2 y_^''"e-[l + i?r/[7^(«,0)]]+e— [l + £^^/[7^(-«,0)r ^ ^^^^ 



where 



7e(u,t) = W^^-^(a + ^^ + zt). (89) 
2 qi- q2 



This result is to be compared with the replica-symmetric result evaluated as 

t 



cosixot) - Erfi{-j=)sin{xot) 



— e{t) (90) 



Both results are compared with each other for short times in Fig, (pi]) and for long times in Fig.(p^). 
For large times we obtain 

GW^^it -> oo) ^ (91) 
with an amplitude smaller than the one of the replica-symmetric solution (|90|), which gives 

G(°)-^(t -> oo) ^ -^./If!!^ (92) 



The behaviour at large times of this T = 0-quantum theory is obviously marked by RPSB. As will become clear below 
the susceptibility x, on which the 0(i)-term of the fermion propagator depends, becomes smaller and smaller with 
increasing steps K of RPSB and finally vanishes in the presumed exact solution at A' = oo. The oscillations similarly 
become slower and disappear at if = oo. 

One may in this context recall the correspondence between long-time behaviour ia, the case of classical Glauber 
dynamics, which was seen to correspond to the nontrivial part of the Parisi functiontHI. 
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FIG. 12. Long-time decay of the fermion Green's function showing gap-induced faster (ORPSB) and slower (IRPSB) os- 
cillations. The amplitude c of the | long time decay is reduced by 1-step RPSB. Oscillations and amplitude c vanish as 
{K cx))-step RPSB (not shown). 



B. Higher order RPSB and 7^ = 00 



The result in 1-step RPSB, although still unstable towards higher RPSB, can be viewed as a much better approxima- 
tion than 0-RPSB, since it already contains features of the full parisi solution at K = oo. In this section we consider 
the RPSB equations at arbitrary K. 



1. K-invariant ratios 



We now derive two K-invariant relations. The first one, between zero temperature gap and nonequilibrium suscepti- 
bility both calculated in K-th order RPSB, is contained in 



(93) 



In order to obtain the second invariant we compare the T 0-limit of the selfconsistent equations for q ~ qi and for 
the DoS at arbitrary order K of RPSB. 

The Green's function at K-th order RPSB can be written as 



kJ''^'Ik^''-L^' kC^^ 



where 



(94) 









'^K-l -1 
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JK-1 


JK-2 Jl 





while Jj^ and C have been defined before by Eq.(65) and Eq.([75|) respectively. 
In the zero temperature limit Eq.(M) turns into 



Gie) - / Ck+i\ 

IK+l 



C-ft:|r=o 



K 



(95) 



(96) 
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The density of states at the gap edge Eg^'^ is hence given by 

,T^^ f f f°° du e 2(92-93) g 2(91-92) 

p(^)(i?f ) + 0) = / a+i|..„ / ai..o - / ^ . . (97) 

00 



A similar expression is obained for the T = gap energy Eg^'^ , since 

r duae'W^ hm^^„/3/^^duie-^I^r^cos/i"'"'(/3ui) 

^ J — oo ^ 



Thus we obtain the second K-invariant ratio 



2 ^^^^ 

Together with the fact that, for each given K, the gapwidth equals twice the fermionic nonequihbrium susceptibihty 
X, which turns into x — PiQ ~ 9(1)) differing only by exponentially small terms from the spin-space-analog x — 
— q{l)) ~ Tu, where denotes the Parisi function q{x) at x = 1, the DoS-hardgaps at finite K terminate in 
a softgap for K — > oo. Note that we did not have to evaluate the T = 0-Parisi function q{x) in order to reach this 
conclusion. Assuming that the relation between gapwidth and x remains valid (at least in good approximation) for 
finite-range models, fluctuation effects should harden the gap. This requires further analysis. 
While we have proved the existence of a spin-glass hardgap at any finite K > with 

Sp^^\T^O,E)^\E-E^g^^\ (100) 

for \E\ > -Eg^' and K finite, the pseudogap of the infinite-range model leads to 

p{T^O,E)r^\E\'' (K^oo) (101) 

with a scaling exponent x, which could eventually become different from one and remains to be determined. This 
pseudogap together with.x=l would be slightly reminiscent of the exponent found for a superconducting glass unitary 
nonlinear sigma modelEj. However exponents differing from one aije-|known even in mean fields solutions for gapless 
superconductors and also at our spin glass thermal tricritical pointE2l. We remark that the pseudogap-solution given 
for the fermionic spin glass model refers precisely to p, — whereas the hardgap-solutions at finite K happened to 
be stable within finite intervals |^| < x/2 corresponding to half-filling only at T =t-0. The regime beyond half- 
filling, identified as the domain of phase separation in the replica-symmetric solutiorJl^, requires further analysis at 
T = as well as several, metallic-, Kondo-type-extensions. Given the experimental facts about phase diagrams in 
HighTcSuperconductorscj'ca the allowance for superconductivity in fermionic spin glass models is interesting too. j— . 
In more difficult models than those treated by us here, the new method of Fourier-Transformations in replica spacec3 
is hoped to facilitate further insight into the highly nontrivial K — oo-solutions. |— . 

We finally note that an overlap distribution function for data clustering shown inE3 and interpreted as a pseudo 
T = 0-problem in a classical spin analogy, revealed, apart from the ratio discussed in Eq.(p9[), a remarkable similaxifaz. 
with the {H = 0)-density of states. It appears interesting to explore pseudo-(r = 0) neural network problemsExEEl 
as potential classical partners of fermionic spin glasses. 



VI. WARD IDENTITY 



Conservation of charge and the corresponding continuity equation provide an exact relation which has a bearing on low 
energy two-particle excitations, connecting them also with the one-particle excitations discussed before. Its feature, 
which is to connect (directly or indirectly) one-fermion propagators with vertices and higher order correlations, is 
particularly important in the present context, since RPSB-breaking is thus transferred into many-partice Green's 
functions. One may not be surprised to find RPSB in spin- or charge-overlap correlations like [(^°'{T)a'' {T')]av for 
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example, but the less obvious strong effects derived before for replica-diagonal quantities will also be present. 
The Ward identity for charge conservation will of course be most important for metallic models, since the replica- 
diagonal continuity equation 

d^v{^} + ^^{J2^UmAt)) = (102) 

where div{j^^} means the discretized divergence of the current operator which is given in terms of the fermion fields 
by 



j-(r,) = -ter, ^ - i'l.t,^^^^ (103) 



links current- and charge-correlations. 

It does not seem necessary to invent replica-overlap currents and densities, if not for the purpose of studying the 
existence of replica nondiagonal Green's functions G""^ involving nonhermitian realizations. It is clear that beyond 
replica permutation symmetry on the level of spins there is the possibility to write, formally, another RPSB on the 
level of fermions. This is however not the subject of the present paper. 

Let us recall the two-fermion propagators K°'^{q,LLj) and the density response function, which played a crucial role 
in the replica-symmetric theory of Anderson localization. Hopping disorder related diffusion of electrons was closely 
linked to the Ward identity. Even if we would not question the stability of those theories against RPSB or perhaps 
other kinds of RSB, the theory of localization due to a random interaction will as well rely on the Ward identity. 
The Ward identity for the insulating model can be viewed as the k — 0-part of the one for a metallic spin glass. 
Ignoring fermion momenta (anyway absent in the insulating case) one may employ 

iu Ap{k ^0,f. + Lu,e) = g{e + lu) - g{e) (104) 

to get 

lim^^Q uj limk^oAp^{k, e + u),e) = 27ri p(e, {qr — Qr-i}) (105) 

shows that the Parisi form of the DoS, depending on all qr — qr-i or on q{x) for K = oo, also enters the vertex 
function. As usual the the three-point function involving the four-current vertex < rT-[a|^(r)ai/(r')j^v/(0)] > defines 
the two-legged quantity Ajti . 

The Ward identity for metallic or superconducting systems with allowance for spin glass order and RPSB has the 
usual form 

ikAj{p + k,p;€ + Lj,€) + iujAp{p + k,p;e + u;,e) = Q{p + k,€ + uj) — Q{p, e), (106) 

where Q denotes again the disorder averaged and replica-diagonal Green's function, which incorporates even at T = 
strong RP SB-effects. 

It is very tempting to mention the consequences of a representation (^4|), although we do not yet know, whether it 
holds beyond the insulating model. It appears likely that it may hold within the Q-static approximation (since the 
Parisi type average given here involves quantities that are static in the insulating case but become time-dependent 
in the metallic case). Thus, within the limits of its validity, this representation could be generalized in a rather 
simple way to the metallic or superconducting cases at zero temperature, and one would obtain the [T = 0) Ward 
identity in the form similar to the one of clean systems but including an disorder average with constrained (cutoff) 
Parisi type distributions (depending on te number of RSB-steps). Of course all order parameters would have to be 
determined selfconsistently for each particular case. Inferring the form of the K =1 propagator, Eq. (p^) , the vertex 
(with external legs) A depends on the Parisi order parameters as (in shorthand notation) 

■ . / , ■ , N e(go)e°^° Es^±i{,,+i^_ - ,,,, 

«'?MAM(P + fc>P;e + '^>e) = / ,G \~, (107) 

This will also occur in metallic spin glasses, whence diffusive modes and conductivity are expected to depend on Parisi 
symmetry breaking. 
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VII. MESOSCOPIC FLUCTUATIONS IN QUANTUM SPIN GLASS PHASES 



Spin glasses have been discussed in the context of mesoscopkisystems0. Nanostructuring of materials with frustrated 
magnetic interactions promise-a large field of applicationsEj. We contributed a recent example of a semimagnetic 
layered semiconductor systemEll. 

Mesoscopic fluctuations have been studied intensively for electronic systems with random potentials. The possibility of 
RPSB was ignored or assumed to be irrelevant or nonexisting. We provide here the technique that helps to reconsider 
this question. This section however briefly evokes only a few examples of mesoscopic fluctuations in spin glass systems. 
A typical example is given by statistical fluctuations of the density of states. Multi-valley correlations described by 

^ PaAEl)PaAE2).-PaAEk), (108) 

where the denote distinct replicas, do not factorize within the spin glass phase. For the choice \Ek\ = E < Eg{H) 
we obtain 

TDt k jE k jE \ 

(109, 

with = i^jj^ and Beta-function B. These multi- valley correlations of arbitrary order k also vanish exponen- 
tially everywhere in the gap as T — )■ and, much weaker than the averaged DoS though, they increase as either gap 
edge is approached at fixed low temperatures. 



VIII. APPENDIX A: QUANTUM FIELD THEORY AND DECOUPLING PROCEDURE 

A. Al: the fermionic Ising spin glass 

The gaussian average over random fluctuations in the magnetic coupling Jy renders the thermodynamic potential 
— Tre"^^^^^^-* of the fermionic Ising spin glass as 



a<b 



(110) 



where 

—b 

A=±l A' = ±l 



We decouple this eight-fermion correlation by a matrix field Q"'" (r, t'). Since the fermionic Ising model does neither 
have spin- nor charge-quantum dynamics, and since Q represents a coarse-grained spin overlap field with the same 

average as the X-field, its static and spatially homogeneous saddle point solution, which has to be determined 
selfconsistently, is extracted by 

Q^'"'(r, t') = A"^'"^' + (5Q^'"'(r, r') (112) 

where the saddle point matrix Q^p = A turns out to be of the usual Parisi form plus a diagonal matrix taking care of 
(cr°)^-averages, hence 

The infinite-range fermionic Ising spin glass is then described by 

n = linin^o- [< Zn > -1] (114) 
n 

with 
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The saddle point contribution for the infinite system {N oo) can be extracted by 



N 



-(/3J)2^(q-) 





exp{ 







a<b 



a<b 



(115) 



, (116) 



while the fluctuation part is expanded in terms of the fluctuation fields SQ. Time-independence of the spin correlations 
for the Ising case including charge interactions and grand canonical description, allows to perform time integrations 
on the fluctuation fields. Hence only the 5Q{uj — 0,uj — 0) components remain. This is different in the metallic case 
and the separation of a dynamic saddle point together with the dynamic fluctuation theory in terms of fields SQiuj, uj') 
was described int3. 

The second class of decouplings, needed to reduce 4th-order Grassmann field products to the integrable 2nd order 
products, must take into account the RPSB in the Q~matrix. With each step K of RPSB te number of new decoupling 
fields increases. The procedure is clearly seen for example at if = 2, where 



< Z„ > = e" 



(^,e^) + i/32j2Tr(Q X) 

^ Parisi 



(117) 



n 



G 

(°3) 



E 



E E 



Q2 = (ct3-l)^ + l cti = (ct2-l)^ + l cto = (ai-l)^ + l 



where the decoupling fields carry a Parisi block index fc (summations over run over replicas pertaining to this 
block) and =9, = y". The generalization to arbitrary K is now obvious. Generating fermion fields 77,77 can be 
added to allow for the calculation of Green's functions. Let us define 



n(n 



n 



) 



n/ n 

(ofc+l) Qfc ,(°fc) ao ,(=0) 
fc+i ~k ~a 



E 



afc + l - 



E 



E 



afc = (afe+l-l)^^^^|fi + l Qfc-l=(Qfc-l)- 



■ E • 

ao = (Qi-l)^ + l 



(118) 



(119) 



The generating functional for an arbitrary number K of RPSB-steps is then obtained as 



2,1(7?, 7;) = e" 

Exp 



(120) 



K 

n 



E 



where the statistically fluctuating Green's function g is given by 

g-i(Q,K-'}) = [5„-i(Q)+^^({4"^'}r^ 



(121) 



The bare propagator g(ti) = l/(ie; + /i) becomes spatially dependent in metallic models, while all of the Ising spin 
glass complications including Parisi RPSB are taken care of by the fields z^. A shift of the fields leads to the result 



Exp{<i>irj,r])) = (11 / d^i')Exp{Y, 9a\^iAz^Wa + v4' - V^) 
= i?xp($(0,0))i?xp(-^7?:'5.(ei,{z^})7;^'). 



(122) 
(123) 



which shows that the averaged Green's function is obtained by the given K+1 gaussian integrations over the fluctuation 
fields, contained in the function g, according to 
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(124) 



Apart from the gaussian weight indicated by the upper index G, one needs the regularized (provided the continuous 
time forniahsm is used; an example for the application of the discrete time formalism is given in Appendix 2 for the 
local Hubbard limit) solution of 



Exp MO, 0)] = l[l[[ef + ^i^ + H^izo,Zl, 



ZK, 



{fcc} 



Y\ \2cosh{l3H{zQ, zi, zk)) + 2cosh{f3^i) 



{bo,} 



(125) 
(126) 



B. Appendix 2: The local Hubbard limit with random U 



The mapping between the Green's function of the fermionic Ising spin glass and the one for a Hubbard interaction 
with random positive U and properly chosen mean value of U, both taken at T = and at half-filling is surprising in 
several respects. A search for more general relations is necessary. A general method is of course the comparison of 
the field theories. A closer look on the field theoretic representation of a random U Hubbard interaction with cutoff 
gaussian SU distribution is needed. 

The decompling of the Hubbard interaction, using the operator identity n^ni — j{{^\n\)'^ — (cr^)^), was given by 
VoUhardtlia. Let us first reconsider the Grassmann integrals using the discrete times (introduced by time slicing) , since 
the Hubbard interaction is one of the nice examples that allows to carry through a regularisation-free formulation. 
The partition function reads 



M-l 



CT fc=0 



A/-2 



exp 



Af-2 



fe=0 
M-2 



(127) 



k=0 a 



k=Q a 



where infinitesimal time steps e = -p- are employed. Using the Grassmann equivalent of the operator identity given 
above, the interaction term is decoupled at any time instant by 



da 



=e sTfce 



/2n 



(128) 



As for continuous timetJll all integrals can be performed exactly. The bilinear exponent is expressed in terms of a 
matrix (S^')^/ ) whence 



Z = liniM^ 
The determinant is found as 

detB{U) 



lim 



M- 



n 



G rG 



detB 



En+nn 



1 + e/i + \/ — {ioLk + A7fc) 



(129) 



(130) 



The remaining integrations over the charge- and spin-decoupling fields are readily evaluated to reproduce the result 
Z = 1 + 26^^^ + e'3(2A'-^) for each site. 

This result being trivially derivable within the operator formalism, the filed theory obviously complicates the deriva- 
tion. Nevertheless field theories provide the best way to find general relations, for which we are looking here. In fact, 
since we know that the infinite-range fermionic spin glass cannot be represented at all temperatures by one zero- 
dimensional random U Hubbard model, we are looking for the differences in those models field theoretic descriptions. 
Before discussing the effect of cutoff gaussian distributions in replicated extensions of the given partition function, let 
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us the generating functional in the discrete time formahsm. 
The generating functional 



^(77,77) = h 



G rG 



-'S'B'S'+riili — 'q-4i 



= Ir 



Y[ detBt 



only requires to invert matrix B. This yields 

Q{r, n) = -^l[r Til B[%,detB- 

^ I. Jak J-lk i^u. 



This is readily evaluated to be (for Tf > ti) 

Q]{Tf-Ti)^--lmM-^oo'W Yl 

I n •J at •J Ik 7. 



k=0 •">''= fci 



1 + 6/1+ V-^(*"fe+7fcJ 



/Z 



(T.k 



M-1 



(131) 



(132) 



(133) 



Since we have now at hand the field theory which correctly reproduces the known results for Z and Green's function, we 
know that the replica formalism applied to a random U local Hubbard limit with the incomplete gaussian distribution 
obtained in ( p^ ) must reproduce the average of In Z and G. It will be interesting and necessary to study in more 
detail the field theory of the random Hubbard model with an incomplete gaussian distribution and its relation with 
the metallic spin glass. 
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